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Recent experiments on unzipping of RNA helix-loop structures by force have shown that ~ 40- 
base molecules can undergo kinetic transitions between two well-defined 'open' and 'closed' states, 
on a timescale ~ 1 sec [Liphardt et al, Science 297, 733-737 (2001)]. Using a simple dynamical 
model, we show that these phenomena result from the slow kinetics of crossing large free energy 
barriers which separate the open and closed conformations. The dependence of barriers on sequence 
along the helix, and on the size of the loop(s) is analyzed. Some DNAs and RNAs sequences that 
could show dynamics on different time scales, or three(or more)-state unzipping, are proposed. 
PACS: 87.15.-v, Biomolecules: structure and physical properties 



I. INTRODUCTION 

Helix-loops are the basic secondary-structure elements of folded single-stranded nucleic acids 
(ssNA). Recent physical studies of single helix-loop RNAs have revealed that despite their simple 
structures, they can display interesting dynamics |l|. When placed under moderate tensions « 
15 pN, telegraph-noise-like 'switching' behavior can be observed. The characteristic time of this 
switching has been seen to be on the ~ 1 sec timescale, surprisingly large given the small size (~ 10 
nm) of the molecules. The purpose of this paper is to present a simple theory capable of reproducing 
these slow switching kinetics. 

The model we use for the energy of a helix-loop structure considers states of partial 'unzipping' 
Q. It is based on available quantitative descriptions of base-pairing interactions and single-strand 
nucleic acid elastic response. In general a large free energy barrier between open and closed states 
for helix-loop structures is found, which qualitatively explains the observed two-state switching 
kinetics. We then combine this model with Eyring-Kramers transition-state theory to access the 
kinetics, and show how this barrier requires a long timescale to be crossed. Finally, we discuss how 
essentially the same approach can be used to describe branched-helix structures, and we predict 
three-state-switching for a specific molecular architecture. 

II. MODEL OF BASE-PAIRING AND SSNA ELASTICITY 

Our description uses the series of molecule configurations obtained by successively breaking base- 
pairs, starting from the molecule ends [^] . In addition to the base-pairing free energy of the double- 
stranded part of the molecule, we take into account the elastic response of the extended, unpaired, 
single-stranded part of the molecule ^[j^. The free energy change associated with opening base 
pair i is simply the difference between the free energies of the paired region, and the free energy 
associated with extension of the open region, 

Ag{i,f)=go{i)~2g,if,t) 

Here go{i) is the free energy of opening base pair i (between 1 ksT and 5 k_gT), obtained using the 
MFOLD program Q . The stretching contribution, gs{f,i), is the free energy per base at constant 
force for stretched ssDNA, which provides a reasonable estimate for ssRNA. Integration of an em- 
pirical fit to the experimental ssDNA force-extension gives: Qs^'^^if) = ksTlss/d ln[sinh(«)/«] 
with u = df/lksT] and parameters Iss — 5.6 A, d = 15 A. Internal unpaired bases (Fig. 1) are 
considered to open along with the base pair immediately preceding them; those base opening steps 
therefore pick up a multiple-base ga contribution. The free energy to unpair the first n base pairs is 
just the sum of the free energies of the individual base-opening steps, 

n 

G{n,f)=J2Ag{i,f) 

i = l 

We emphasize that both the base-pairing and the elasticity contributions to (|l|) and J^) are free 
energies, i.e. are coarse-grained over atomic-scale fluctuations. The above zipper model pTis roughly 
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equivalent to that described in the Supplementary Materials of the presence of the stretching 
force allows to discard complex opening/closing pathways relevant at zero force ^ . Note that some 
cooperativity between base pairs is introduced by MFOLD (the pairing free energy of neighboring 
base pairs is not simply additive, but includes some sequence specific stacking interactions). 



III. EQUILIBRIUM BEHAVIOR OF A SINGLE HELIX-LOOP STRUCTURE 

Fig. 2a shows G{n, f) and the equilibrium probability P(n, /) cx exp(— Gfp, f)/kBT) for a simple 
RNA hairpin, called P5ab (Fig. la), under the solution conditions of [Ql near the critical force 
/* ~ 15 pN where two-state switching is observed. At the critical force, the open and closed states 
(n = 22 and n — 3, note n is just the number of broken base pairs) dominate; below or above 
this force, the free energy landscape is tilted either to favor the n = 3 or n = 22 state. At the 
critical force, there are various barriers due to drops in free energy resulting from openings of the 
U bulge, the weak non Watson-Crick GA central pairs, and the final loop. The largest barrier ~ 11 
ksT must be crossed to reach the closed state from the open one, and vice- versa. The two-state 
behavior observed for P5ab follows from the partition of probability into two peaks separated by an 
improbably-accessed barrier region (see Supplementary Materials, Ref. Q). 



IV. DYNAMICAL MODEL 



To reach a quantitative understanding of switching, we introduce a dynamical model for the 
motion of the 'fork' separating the base-paired and opened regions of the molecule We propose 
the following expressions for the elementary rates of opening and closing base pair n {i.e. to move 
the boundary between the open and closed portion of the molecule from n to n — 1 or to n 1): 

r<,(n) =r e-»«'"''"=«^ , re(/, n) = r e-^''='-f'"'/'=«^ (3) 

Here r is essentially the microscopic rate for a base pair to move together or apart in the absence of 
tension or base-pairing interactions, or roughly the inverse self-diffusion time for a few-nm-diameter 
rrji^/kBT lO'^ s"\ with 1=10 nm, rj = 0.001 Pa sec, and fesT = 4 x 10"^^ J. 
we have made the simplifying approximation that the opening rate Vo has no force de- 
pendence, and is simply proportional to the exponential of the base-pairing free energy of (^). 
Eyring-Kramers transition-state theory applied to breaking of a chemical bond considers indeed 
the potential energy to be 'tilted' by a force-times-displacement contribution. Because hydrogen 
bonds break for relatively small displacements (~ 0.1 nm) the reduction in the potential energy of 
the single-base-opening transition state will be roughly 15 pN xO.l nm = 0.3 ksT. This can be 
neglected with respect to the base-pairing free energy of a few ksT, which is a lower bound to the 
energy of the transition state associated with breaking a single base pair. Detailed balance then 
determines the closing rate Vc to be proportional to the exponential of force times displacement, i.e. 
to the energy of a fluctuation that is able to pull the two bases back together in opposition to the 
applied force. 

The rates ^ lead to a master equation for the probability p„ (t) for the boundary to be at site n 
at time t: 



, , N 
dpn{t) _ _ \ - 

dt ^~ ^ 

m = 



T„,„ p„,{t) (4) 



This {N X (A'^ + 1) matrix Tn,m is tridiagonal, with nonzero entries Tm-i,m = 7'c(/, rn), Tm+i, 

ToijTl).^ and Tm.m — Tjn — \.7n Tm-{-l,m. 



V. SWITCHING KINETICS OF A SINGLE HELIX-LOOP STRUCTURE 

We have solved (diagonalized) (^ for P5ab (Fig. la). The smallest eigenvalue is 0; the eigenvector 
is the equilibrium Boltzmann distribution. At the critical force where the molecule is on average 
half-open, the smallest non zero eigenvalue is Ai = 2.1 x 10~^r, which corresponds to the slowest 
mode of fluctuation, the 'switching' of the boundary of the open region between n « 3 and n « 22. 
The remaining 21 eigenvalues are all well separated from the leading ones. The second largest 
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eigenvalue is A2 = 0.9 x 10~*r (Table). Thus the theoretical dynamics of P5ab involve one slow 
opening-closing transition, combined with many other transitions occurring more than 50 times 
faster. The net rates of opening (fco) and closing (fcc) can be computed from Ai — (ko + fcc), and the 
ratio of the open and closed equilibrium probabilities equal to kc/ko- To compare our theoretical 
results with the experiments of [Q| we have fitted our result for k* = ko — k^ at the force /* where 
the open and closed states have equal probability to experimental data giving r = 3.6 x 10^ 
sec~^. 

Fig. 3a shows time series from Monte Carlo simulations of the molecular motion; slow two-state 
switching is seen on a « 0.25 sec timescale, on top of which occur much faster small fluctuations. 
When we convolve these data with a 20 Hz low-pass filter (as used experimentally the result 
(Fig. 3b) is essentially the same as the experiment. The variation of the rates with force given by 
theory are also in good agreement with experiment(Fig. 4). The 'transition state' coordinate can 
be inferred from the relative slopes of the opening and closing rates of Fig. 4 (independent of the 
fitted parameter r) around the critical force Q. The critical-force transition state is at the top of 
the free energy barrier, at n* =11 — 12 (Fig. 2a). 



VI. BARRIERS FROM GENERIC HELIX-LOOP STRUCTURE 

For molecules with repeated sequences and no terminal loop e.g. AU followed by a long GC helix or 
a crosslink between end bases, the kinetics essentially corresponds to a diffusion in a flat free energy 
landscape, and shows none of the two-state character of the P5ab RNA. For a repeated AU sequence 
of 71 = 25 base pairs, the switching time is just the diffusion-like time t* ^ 2iVV(7rV*) = 2 X 10~* 
sec (Table). Thus, the IQQQ-times longer switching time of P5ab comes from the 11 ksT barrier of 
Fig. 2a. 

The presence of a loop is sufficient to generate such a large barrier and, consequently, two-state 
switching. The simplest illustration is given by a homogenous DNA sequence ending with a loop. 
We have considered a 24-base Poly(GC) homogeneous-sequence helix terminated with a 4-base loop 
(Fig. lb, 2b) Our theory predicts a two-state switching behavior on the same time scale as 

P5ab. The switching time is 50 times larger for a longer 8-base loop (Table). 

The free energy barrier G* at criticality for a S-base-pair stem (uniform pairing free energy go) 
followed by a L-base loop (closing free energy gioop{L) at zero force) (Fig. 1) can be simply estimated. 
The critical force /* is given by the condition that the free energy of the open molecule equals the 
free energy of the closed molecule, G(0, /*) = G{S, /*) = S go - (23 + L) gssif) - gioop[L). The 
barrier height G*{S, L) = G{S — 1, /*) then reads 

G*(5,L) = (5- l)(,o - 2,..(r)) = (S-'n9o^Y^l9.o.AL)) 

Table 1 shows that, for a fixed stem length S, the critical force decreases with the length L of 
the loop, while the free energy barrier G*, and the switching time t* increase. For non-repeated 
sequences, (^) is only approximate when substituting go with an average pairing free energy; it 
allows an estimate of how the switching time depends on S and L . Note that the critical barrier 
essentially depends on the smaller of the two lengths S, L. 



VII. BRANCHED-HELIX MOLECULES AND MULTIPLE-STATE-SWITCHING 

Our approach can be extended to more complicated situations e.g. nucleic acids with branched 
structure. An example is the P5abcAA RNA molecule Q of Fig. Ic, with free energy landscape as 
shown in Fig. 2c. The opening of the first 12 bases (helix HI) follows as above, but going past the 
branch, description of the independent opening of the two helix regions requires a three-dimensional 
free energy landscape i. e. free energy as a function of the positions of the two unzipping boundaries 
]l^ . A rich behavior emerges, shown in Fig. 3c. At the critical force /* = 12.9 pN, HI switches 
on a long time scale t* = 10 sec (Table), while the short lateral helix (H2) opens and closes with a 
much shorter characteristic time t2 — 9 msec. The lateral long helix (H3) opens very rarely at this 
force. Predictions for the rates, barrier height, ... are reported in Table 1. 

It is possible to design molecules with multiple-state dynamics. Consider the molecule of Fig. 
Id, which has two well-bound Poly(GC) regions separated by an unpaired bubble, and terminated 
by a loop. For this system, three-state switching occurs (Fig. 3d), and should be observable in the 
frequency range accessible in experiments such as that of 111 . 
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VIII. CONCLUSION 



A few improvements might be added to this model. First, opening may occur through the 
cooperative nucleation of a few base pairs bubble [Q. Second, mismatches might take place during 
closing ||l^ , although they are highly limited by the presence of the 15 pN force. Finally it would be 
interesting to be able to describe unzipping events involving breaking of tertiary structures; these 
are thought to be present in the molecule P5abcAA in presence of Mg^"^ jllJl^ . 

A general result of our work is that slow switching character should be quite generic for small 
biological molecules with helix loop structure. We note that our approach could also be used to 
analyze the opening-closing dynamics of nucleic-acid-detecting DNA 'beacons' [^, both on their 
own, and in the presence of their targets. In this case the 'unzipping' forces are applied by the 
hybridization interactions instead of by a large force transducer. Since such experiments amount to 
molecule recognition processes it is not impossible that slow barrier-crossing transitions of the sort 
discussed here occur in-vivo. 
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TABLE 1. Theoretical results for the molecules of Fig. 1 compared to experimental values from [1, 14] (in bold) when 
available. Columns indicate the free energy AG — G{N, / = 0) at zero force, the critical force /*, the rate of opening/closing 
k* , the switching time t* and the free energy G* of the highest barrier at criticality, the variation of the opening and closing 
rates upon force, and the ratio of the two largest non zero eigenvalues. This ratio is small when open and closed states are well 
defined, and close to one otherwise. Uncertainties in base-pairing free energy are at most Sg ~ 0.5fcflr/bp, with consequent 
total uncertainties of « N^^'^Sg (or about SksT for N = 25) for AG and of ^ [N/2Y''^5g ~ 2fcsr for G*, and « 3fcsr/20 
nm« 0.6 pN for the critical force. Notes: " the ratio is close to one due to the presence of three, and not two, states giving rise 
to two large barrier crossing times (the remaining fluctuation times are much shorter: A1/A3 ~ 0.0001); " the predicted rate 
(0.1/sec) for P5abcAA is very close to the lowest frequency (0.05 Hz) resolved experimentally [1], thus its value is known only 
roughly. 
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FIG. 1. Nucleic acid unzipping experiment and molecules studied. A constant force is applied to a helix-loop structure while 

the distance between molecule ends is measured, (a) P5ab RNA, a single helix-loop structure present on the P4-P6 domain of 
a self splicing group I intron of the Tetrahymena thermophila. The structure shown is predicted by Mfold [7] apart from the 
G-A weak pairs [10] indicated here with dots, and with the U-bulge translocated; (b) DNA hairpin consisting of a poly(GC) 
helix terminated with an ATAT loop [14]; (c) P5abcAA RNA, a variant of P5ab with an additional helix giving a Y-branched 
structure at zero force; (d) Hypothetical RNA molecule obtained by ligation of two poly(GC) helices as in (b) and replacing A 
bases with Ts. 
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FIG. 2. Free energy landscapes (in ksT, top) and probability distributions for the position of the opening fork (bottom) at 

the critical force for molecules of Fig. 1. (a) P5ab; (b) Poly(GC) RNA with: a 8bp loop (diamond), a 4 bp loop (squares); no 
loop (circles), (c) P5abcAA; (left: opening positions ni = 0, . . . , 12 corresponding to helix HI, right: opening positions (n2, ns) 
with n2 = 0, 5 - H2 - and ns = 0, 9 - H3 -); by definition configurations ni = 12 and (n2 = 0, ns = 0) coincide. 
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FIG. 3. Unzipping kinetics at the critical force. Distances between extremities of the molecule are shown for: (a) P5ab; 

there is a slow switching between the n ~ 3 and n ~ 22 configurations and fast transitions between configurations n around 
these ones (Fig 2a). (b) Convolution of (a) where oscillations faster than 20H[z are averaged out; (c) P5abcAA; there is a slow 
switching between the closed molecule and the molecule with HI opened and H2 opening or closing on a shorter time scale, 
the opening of H3 (distance between extremities of 31 nm) is a rare event at the force of 12.9 pN; (d) the hypothetical RNA 
molecule of Fig. Id. 




FIG. 4. Log. of the opening and closing rates for P5ab as measured in [1] as a function of force (full circles: opening, 
empty circles: closing), compared to theory (full line: opening, dashed line: closing). The slopes of In fco. In fee give the relative 
positions rio = 8, nc = 11 of the transition state from the closed (n = 3) and open (n = 24) states respectively, with an absolute 
location in n* ~ 12. 
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